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Abstract. Log-periodic amplitudes appear in the critical behavior of a large class of 
systems, in particular when a discrete scale invariance is present. Here we show how to 
compute these critical amplitudes perturbatively when they originate from a renormal- 
ization map which is close to a monomial. In this case, the log-periodic amplitudes of 
the subdominant corrections to the leading critical behavior can also be calculated. 



Dedicated to Herbert Spohn, our colleague, friend and source of inspiration, 
on the occasion of his 65 th birthday 

1. Introduction 

Power laws with non integer exponents are common in the study of scale invariant struc- 
tures, in particular in the theory of critical phenomena at second order phase transitions. 
In some cases, when the scale invariance is discrete [HI [361 CEH EH USE IS [21 Q], the 
amplitude of the powerlaw acquires a periodicity, often called log-periodic oscillations: see 
[351 [21] for reviews on the topic. These oscillatory amplitudes are usually more difficult 
to calculate [131 \12\ El [26] than the critical exponents. When the scale invariance is 
associated to a renormalization group map, they can be related to some properties of the 
map, for example to the shape of its Julia set, see [131 IB]. 

The goal of the present paper is to show that these oscillatory amplitudes can be com- 
puted perturbatively when the unperturbed map is simple enough, in our case a monomial 
map. To be more precise let us consider the following question: take the map x h-» F(x) 
with 

F(x) = F e (x) = x p (1 + e G(x)) , (1.1) 

where p = 2, 3, . . . and the perturbation G{x) is sufficiently regular and has a finite limit 
for x — > oo (we will actually need a condition on the decay of G'(x) too, and this will 
be pointed out later on, see (14. 7p ). If x* is the largest real fixed point of the map (so 
that )), this fixed point is unstable, that is F'(x*) > 1, for e > sufficiently 

small. We will actually just focus on F(x) for x > x*. If F^ n \x) denotes the n-th iterate 
of the map F, the sequence ip n {x) defined by 



logi?K>( 



Ms) = „„ , 

has a limit for all x > x* 



P" 



ipoo(x) = hm ip n (x) . (1.3) 
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One directly sees that i^ooix) > for x > x*, while ipoo(x*) = 0. The existence of the 
limit directly follows from 



\^ n+1 {x)-^ n {x)\ <p- n -'log 1 + esup |G(x)| . (1.4) 

V x>x* J 

It is straightforward to see that the limit satisfies the functional equation 

p^oo(x) = i/j oc (F(x)), (1.5) 

which goes under the name of Bottcher equation in the iterated function literature j3[ [30] , 
and one expects and sometime can prove - see for example [SJ [H] - that 

^w^f-^K-^f)- (L6) 

where the exponent a is usually non-integer and B(-) is a strictly positive periodic function 
of period one. The exponent a has the well known expression 

a logF'(x*)' 1 ' 

while the oscillatory amplitude B{y) is usually hard to determine. One goal of the present 
paper is to show how to compute this oscillatory amplitude perturbatively (related small 
e expansions can be found in [9]). We will see in Section [3] that all the subdominant 
corrections to (jl.6p can be calculated as well since they all can be expressed in terms of 
the periodic function B(-). 

Remark 1.1. (jl.6p is easily guessed: if one assumes that ^oo(x) ~ (x — x if ) a C{x — x*) as 
x \ x*, with y i-» C{y) defined for y 6 (0, j/o) for an arbitrary choice of yo > 0, continuous 
and such that < C m \ a < C{y) < C max , one gets from (jl.5p that 

(x-x,) a C(x-x„) ~ (x - x*) a c(>(x,)(x - x,)) , (1.8) 

and this implies the value (|1.7|) of a and the fact C^F' (x^y^j = C(y) for every y, 

that is that C(-) is multiplicatively periodic of period F'(x^). Equivalently, C(y) = 
B(—log(x — x*)/logF'(x*)), with B(-) periodic of period one, as in (|1.6j) . 

Remark 1.2. In what follows we will often assume that 

G(1)=0, (1.9) 

so that x* = 1 when e is sufficiently small. There is no loss of generality in doing so 
because the general case can be reduced to G(l) = by the change of variable x i— > x/x* 
and easily be computed perturbatively in e 

^.-.M^ ^yw ,,. (1 . 10) 

p — l z(p — iy 

In this paper we will also discuss the case where the map F(x) has the form 

F(x) = F p (x) = qx p + (1 - q)P(x) , (1.11) 

where p is a large positive integer, q € (0, 1) and P(-) is a polynomial such that -P(l) = 1. 
We aim at understanding what the oscillatory behavior becomes for large p (when P(-) and 
q are fixed). The assumptions guarantee that, for p sufficiently large, 1 is the largest fixed 
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point (and it is unstable). Maps in this class are for example (x p +l)/2 and (x p +2x 2 — x)/3. 
We will see that for large p the oscillatory amplitude can be sharply estimated. In a 
nutshell, the reason why the map becomes amenable for large p is that sharp estimates 
of the iterated map are easy as long as x is very close to x* - the linear regime 

- and when x is large, because the map there is close to qx p . The only difficult part of 
the iteration is in between these two regimes: we will see however that for p large this 
intermediate regime essentially reduces to one iteration which one performs explicitly. 

2. TWO EXAMPLES 

Let us discuss here briefly two different examples where the above question is relevant. 

2.1. The wetting problem on a hierarchical lattice. One of the developments in 
the theory of critical phenomena by the renormalization group approach was to study 
statistical mechanical models on hierarchical lattices [H [28| [231 f!5| [31] . These lattices 
are usually constructed by a recursive procedure which is at the origin of a discrete scale 
invariance and which allows one to write exact renormalization transformations. 

Here we consider the so-called diamond lattice |28[ [15] (in fact a generalization of 
it) which is constructed as follows. One starts with a lattice composed of a single bond 
between two sites A and B. Then at every iteration step each bond of the previous step is 
replaced by p x b bonds, namely by b parallel paths of p bonds each. Therefore if, at step 
0, one starts with a single bond between A and B, then, after n iterations, there are (pb) n 
bonds in the lattice, 6^ pn_1 ^^ p_1 ^ directed paths from A to B, and each of these directed 
paths has a length L n = p n . 

One of the simplest critical phenomenon one may consider on such a hierarchical lattice 
is the wetting transition |11 [ I19 [ [20]: in the wetting problem, all the bonds along a single 
path (the special path) from A to B have an energy —a whereas all the remaining bonds 
on the lattice have energy 0. The energy Eyy of a path W from A to B is simply the 
number of bonds it has in common with the special path. Equivalently this energy E\y is 
the sum of the energies of the bonds visited by the path W. Then the partition function, 
at temperature T, is defined as 

Z ( T) = £exp(-^V (2.1) 

where the sum is over all the b p- 1 paths W. 

When the temperature varies, one observes a transition between a bound phase, where 
typical paths have a non zero fraction of their length in common with the special path, to 
a free phase where this fraction vanishes. At infinite temperature, all paths have weight 1 
and the partition function counts simply the total number Z n (oo) of paths 

Z n (oo) = bl^ = b (Z n „!(oo)) p , (2.2) 

while at finite temperature, the partition function satisfies the following recursion 

Z n (T) = Z n ^{Tf + (b - 1) Z„_!(oo) p . (2.3) 

If one introduces the ratio x n — Z n (T^ j Z n {po) , it follows that 

x p _|_ 5 _ i 

x n = F(x n _i) where F(x) = , (2.4) 
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with xq = exp(a/T). At the n-th step the free energy f n per unit length is 

UM != _ i^i^ + . (2 . 5) 

which implies for the large n limit of /„ (recall (jl.2p and (|1,3I0 



log 6 

/oo(^o) = V ; oo(2;o)H 7- (2.6) 

p — 1 

Note that, while from the construction of the lattice b = 2, 3, ... is an integer, the iteration 
(|2.4p makes sense for every b > 1. In what follows we will consider this generalization (see 
Remark 12.31 for another probabilistic interpretation of Z n in this generalized context). 
Finding the critical behavior of the free energy when the temperature is close to the 
wetting transition temperature, i.e. when xq is close to the fixed point x* of the map F, 
is a problem of the same nature as the one posed in the introduction and the question we 
address here is to try to determine the periodic function B(-) in front of the power law 
singularity 

/.(«) - ^ •> (x - B (-^^f) ■ (2.7) 
p- 1 V iogF'{x*) ) 

The map (I2.4p can be reduced to the form (II. ip or (jl.lip in the two following cases 

i 

• when b is close to 1 or when b is large, by making the change of variable x = bp- 1 X, 
the map (|2.4p becomes 

X -> X p + 6"^ _ &-^t ? (2.8) 

which is of the form (|1.1D . 

• when b is arbitrary and p is large, because the application (|2.4j) coincides with 
([LIB , with the choice q = 1/6 and P(-) = 1. 

Remark 2.1. It is worth pointing out that the partition function of the wetting model 
on more general hierarchical lattices is obtained by iterating a polynomial of the form 
F(x) = qo + q\x + • • • + q p x p , with the qj's probability weights. Of course F(l) = 1 and 
it is easy to see that there is at most another positive root: recall that the unstable fixed 
point a;* on which we focus is the largest positive root. We record here the well known 
formula (e.g. [2U [Til E]) that actually holds for every x € (0, oo) 

, oo 

^oo(x) = l -^+logx + J2p Hl+1) Q(F {i) (x)) , (2.9) 

with Q(x) := log(F \x)x~ p / q p ) . It is rather straightforward to extract from this formula 
that ^oo(') is (real) analytic on (x*,oo) (e.g. [8]): in statistical mechanics terms this 
says that is the only critical point (note that with the definition (jl.2l 1 1 . 3|) one has 

4>oo(x) = for x e (0, x*]). 

2.2. The Galton Watson process. The Galton process is a simple model for the size 
of the lineage of an individual [25] . The generations do not overlap and all the population 
is replaced at every generation. 

One starts with a single individual at generation g = 0. Then each individual living 
at generation g has a probability qt of having k offspring at the next generation g + 1. 
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Therefore if N g is the number of individuals at generation g (so Nq = 1), one has 

N g+l = Y, k i, (2.10) 



where the fcj's are N g independent random variables which take the value k with a prob- 
ability qk- The size N g of the population at generation g is therefore a random variable. 
The average number of offsprings at generation g is (see (|2.10p ) given by 

w = (^ fc9fe ) = (A;)9 - (2 ' n) 

So the expected size of the population grows exponentially for (k) > 1. It is well known, 
see e.g. [25], that if (k) < 1 then after sufficiently many generations the population goes 
extinct. We will actually restrict to the supercritical case (k) > 1 and in this context it is 
natural to define [24^ [25] the generating function of the ratio N g /{N g ) 

H g (X) := ^exp(Ai)^ , (2.12) 

One can see from (|2.10j) that 

*.w - E (w))) k » " F (*-' (w)) = F( " (*' ' (2 ' 13) 



k 



where F(x) = ^2 k qk % ■ Then H^X) := lim 9 _s.oo H g (X) satisfies 

^oo(A) =f(h x (±\\ , (2.14) 

This functional equation is known under the name of Poincare equation [30] and Hoo(") is 
the unique solution of (|2.14p which is analytic in a neighborhood of the origin and such 
that Foo(0) = 1 and fl^(0) = 1 [3D]. Note that that, given #^(0) and ^(0), one can 
extract recursively from (|2.14p all the coefficients of the Taylor expansion at A = 0. 



Remark 2.2. To be precise, the limit exists in great generality and it is (obviously) finite 
if A < |25[ Ch. I]. In fact the random variable W g := N g /(N g ) converges almost surely 
to a limit random variable Woo as soon as (k) < oo: this is particularly elementary if 
(k 2 ) < oo, because in this case one can explicitly compute ((W g + n — W g ) 2 ) for n € N [251 
p. 13] and establish the convergence (in probability) of W g to Woo and the fact that Woo 
is not identically zero. Of course in this generality i?oo(A) can be oo for A > 0, but we 
will focus on the case in which only a finite number of qkS are non zero and it is therefore 
rather straightforward to see that for every A > there exists C\ such that 

sup (exp(AW 3 )) < C x . (2.15) 

9 

This not only implies that ifoo(A) is finite for every A G 8 and that it is equal to 
(exp(AWoo)), but also that the same is true for any A complex: as a matter of fact -ffoo(') 
is an entire function, that is analytic on the whole of C. 
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Let us therefore assume that only a finite number of q^s are non zero, so that the map 
F(x) is polynomial of degree p = max{/c : q^ ^ 0}. T. E. Harris has proven [23] that for 
large A 

with a = logp/ log(/c) as in (|1.7|) and L(-) a continuous positive periodic function of period 
1. Harris was unable to establish that L(-) is not a constant as soon as q p < 1 (for q p = 1, 
that is for -F(x) = x p , it is straightforward to see that L(-) is constant) and this open issue 
has drawn the attention in the mathematical community (see for example E] , dealing 
precisely with the problem left open by Harris) and similar - sometimes strictly related 
- oscillatory behaviors have been pointed out repeatedly (e.g. [271 El E31 E21 [37] and 
[8] for further references: the stress is often on the nearly constant and nearly sinusoidal 
character of these oscillations). To our knowledge, establishing in general that L(-) is 
non constant is still an open problem for p > 2 (for p = 2 the non triviality of L(-) is 
established in [8] by exploiting results in [10]). 

Here again one can try to perturb around a situation where the map becomes simple. 
However there is the following important preliminary issue: what is the relation, if any, 
between the periodic functions B(-) and L(-), both of period one, appearing in (|2.7p and 
in (|2.16p , A priori, these two periodic function come out of different questions that have 
in common only the underlying polynomial F. But in fact if the polynomial is the same 
in the two questions, then the periodic functions coincide up to the change of variable 
B(y) = L(—y) (this has been proven in [8] and the need for switching the sign is due to 
the way we have defined B(-) and L(-): the argument of proof is reproduced below). Given 
this observation, the analysis of B(-) and L(-) in the two perturbative limits we consider 
becomes just one problem. 

2.3. Oscillations for hierarchical wetting and Galton- Watson models. We explain 
now how, by exploiting the two functional equations f| 1 . 5 1) and (|2.14j) . one can see that 
the periodic function B(-) in (|2.7p and L(-) in (|2.16p coincide up to a sign change in the 
argument. For this we fix a polynomial F(x) = X^f=o 1^ > n ^ e m Remark [2J] for the 
wetting model and in Remark l2.2l for the Galton- Watson process, and we define for every 
A > 

M(A) = A->oo (#oo(A)) , (2.17) 
and by applying first (|2.14p together with a = logp/ log(fc), and then (11.5P we have 

M(X(k)) = i^^ oo (F(F 0C (A))) = X~ a i>oo (-Hoo(A)) = M(A), (2.18) 
p 

that is, M(-) is multiplicatively periodic of period (k). Moreover M(-) is analytic on (0, oo) 
because -Hoo(') is entire (cf. Remark [2. 2 p and ?/>oo( - ) is analytic on (1, oo) (cf. Remark l2.ip . 
The fact that M(-) > on (0, oo) is just the strict positivity of ipoo( x ) f° r x > 1 and the 
fact that H 00 (X) > 1 for A > 0. 

Now from the definitions it is immediate to see that ipoo( x ) ~ logx for x — > oo and that 
lim A _>oo floo(A) = oo, so (I2.17P 

log(#oo(A)) A ~°°A Q M(A), (2.19) 

which provides a proof of (|2.16p and relates explicitly M(-) and L(-). 
On the other hand (I2.17P can be rewritten as 

^oo(^) = (H^{x)) a M(H^\x)) . (2.20) 
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But since iToo(O) = 1 and H'^O) = 1, we see that H^il + y) ~ y for y \ 0, so that from 
(I2.20P one directly obtains for x \, 1 

^ooN ~ (x- l) Q M(x- 1) , (2.21) 

which proves (|2.7|) and explicitly relates B(-) to M(-), and we see that B(y) = L(—y) for 
every y. This completes the proof of the equivalence of the two example models we have 
proposed, as far as oscillations are concerned. 



Remark 2.3. As it was pointed out in |2Uj . the partition function of the hierarchical 
wetting model has a Galton- Watson representation: 



2.(00) - \"*(t""))- (2 ' 22) 

Formula (|2.22p shows that the partition function of the hierarchical wetting model can 
be interpreted as the generating function of the size N n of the population at the ra-th 
generation of a Galton- Watson processes. It also provides an expression for the partition 
function of the hierarchical model in terms of sum of Boltzmann weights when b is not 
integer. 

3. Main results 



3.1. Oscillations for the map fll.lHl.9j) for e small. Here instead of writing the results 
in terms of the difference x — x*, we found slightly more convenient to use \og{x/x if ). One 
can of course make a simple change of variables to rewrite the final result in terms of the 
difference x — x*. Then (jl.6p becomes for x close to x* 

log(log[x/x*])' 



^oo(x) ~ (log[x/x,]) a B 



(3.1) 



logF'(x*) 

where £>(•) is a periodic function. By replacing i/joo(x) by this leading behavior in (jl.5 
one can calculate subdominant contributions when x — > x* and one gets at next order 

log (log[x/x*]) 



V'ooOe) 



{\og[x/ Xif }) a B 
- (log[x/x*]) 



log F'(x*) 
+1 a[x*F"(x*) + F'(x* 



F'(x*) 2 } 



B 



2F'(x*)(F'(x* 
log (log[x/x*])' 



-1) 

1 



(3.2) 



-B' 



log(log[x/x*]) 



+ 



\ogF\x*) J logp \ logF'(x*) 

In fact all the subdominant contributions can be expressed in terms of derivatives of 
the map F(x) evaluated at x* and of derivatives of the periodic function B(y). So the 
problem of calculating the amplitudes of the dominant contribution and of all subdominant 
contributions reduces to the calculation of the periodic function B{y). 

One of our main results, derived below under the condition (ll.9p is that 

B(y) = 1 + eJ3i(y) + e 2 B 2 {y) + ... (3.3) 

with 

-i 



Bx(y) 



1 

p 



■yG'(l) + ^G(e 

k=0 



k-y 



+ 



y~kQ 



n k-y 



G'(l) 



(3.4) 
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and 



+ 



E 

fc>0 



-^G 3 ( e P k -y\ +e P k -y G > ( e P k -y 



B 1 {y)-P y - k u(e p 



Q k-y 



E 

k=— oo 



p 2 



B\(y) — P y ~ U ( e 



where is denned by 



1 oo 

U(x) = ~Y,P~ k G(x pk 

P k=0 



(3.5) 



(3.6) 



Remark 3.1. In the particular case 

F(x) = (l-e)x p + e 
one has G(x) = x~ v — 1 and therefore 



Bi{ y ) = y + Y,g(y- k )+ Yl (a(y-k) + i) 



k=l 



k=—oo 



with 



(3.7) 



(3.8) 



g(y) = f (exp(-p-») - l) . (3.9) 

From this expression one can obtain the Fourier coefficients a n = Jq exp(— 2-nniy)Bi{y)dy 
of the periodic function JBi(-) and one finds 

1 l + r'(l) 1 1-7 



2 logp 2 logp ' 

(7 = 0.5772 ... is the Euler-Mascheroni constant) and for 

1 „ / 2iirn 

r 1 



(3.10) 



(3.11) 



log p \ log p 

These espressions are compared in Table 1 with the Fourier coefficients calculated directly 
for a finite value of e. 



Remark 3.2. (|3.10p and (|3.1ip follow from (|3.8p by direct computation. It is however 
worthwhile to point out that there is a direct link with a Poisson summation formula (and, 
for example, with a similar computation in a different context [351 § 3.3]). In fact, we have 

00 

B[(y) = 1+ £ g'(y + k), (3.12) 

A:= — OO 
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1/e 


(oo - l)/e 


K(ai/e)xl0 8 


3(ai/e)xl0 8 


5R(a 2 /e)xl0 14 


3(a 2 /e)xl0 14 


10 


0.129 


-3.65 


-7.08 


1.47 


-0.80 


10 2 


-0.093 


-4.80 


-7.07 


1.62 


-1.11 


10 3 


-0.108 


-4.788 


-7.164 


1.682 


-1.085 


10 4 


-0.0978 


-4.786 


-7.174 


1.6887 


-1.0819 


10 5 


-0.10993 


-4.7858 


-7.1753 


1.68938 


-1.08166 


oo 


-0.10995 


-4.7858 


-7.1757 


1.689459 


-1.081618 



TABLE 1. For F(x) = (1 - e)x 2 + e the limit e \ expressions in (|3.10[1 and in (|3.11ll 
are compared with values for e > obtained by using formula (|2.17[) : an accurate 
expression for -Hoo( - ) is obtained by writing enough terms in the Taylor expansion at 
zero (these terms are easily obtained recursively by using Hao(0) = 1, #io(0) = 1 and by 
taking derivatives of the Poincare equation (|2.14[) ) and by exploiting the series in (|2.9[) 
to get a good approximation of ipoo(-)- 



so the n-th Fourier coefficient of B[( ) is equal to J"_f° g\y) exp(— 2mny)dy =: g'(2irn) for 
(Poisson summation formula, see e.g. |344 Ch. 4]), while the 0-th Fourier coefficient 
is of course 0. It is then straightforward to see that a n = ig'(2irn)/2irn and (|3.1ip can be 
recovered this way too. 

3.2. Large p limit of the map (II. lip . For the map (Jl.lip . for q and P(-) fixed and p 
large, so that in particular F'(l) = qp + (1 — q)P'(l) ~ qp, our result is that if 



log(x - 1) 



log(F'(l)) 

where n{x) is an integer, then for y ^ 1 as p — > oo 



n(x) + y with r\ < y < 7] + 1 



(x - l) a q y for rj < y < 1 , 
(x - l) a q l ~ y for 1 < y < 1 + r] . 



(3.13) 



(3.14) 



This gives a discontinuous amplitude at y = 1. Recalling (|1 .6j) , an equivalent way to state 
(|3.14p is to say that for every y Z 



lim B(y) = q 

p— >oo 



-{y} 



(3.15) 



where {x} is the fractionary part of x, that is {x} = x — max{n 6 Z : n < x}. 



One can resolve the discontinuity at y = by analyzing, in the large p limit, the range 
y — 1 ~ (logp)" 1 and one gets 



^(x) ~ (x- l)° e -( 1 -2') lo gf log 



D (l-y) logp 



(/ cxp 



+ 1-9 



(3.16) 



Note also that in the case q = 1 — e and P(-) = 1 for large p one gets from (|3.4p - (|3.5p in 
the range < y < 1 

1/(1/ + 1) 



#i(y) = 2/ and £ 2 (y) 



(3.17) 



which agrees with (|3.14|) when one takes the first two terms in the expansion in powers of 
c. 
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y 

Figure 1. With the choice F(x) = (x p + l)/2 we have plotted B(-) for 
p = 10 3 (dot-dashed line), 10 5 (dashed line) and 10 7 (solid line), and the 
this dotted line for p = oo. For the p = oo case we have exploited the 
result (|3.15p . that is we have plotted y i— >• 2^ y \ The numerical data for 
p < oo have been obtained exploiting (12 .9p . The convergence is rather slow 
simply because the effective large parameter is log p. Notice that already at 
p = 1000 B(-) is far from being nearly constant and starting from p = 10 5 
it is also evident that it is not nearly sinusoidal. These large amplitudes 
contrast with the small amplitudes one usually observes for small p [5j El 

mug. 



4. Perturbation in powers of e far from the critical point 

Our approach consists first in calculating in this section ipoo(x) in powers of e, for x ^ x*. 
Then by matching (in Section [6} the singularities (obtained in Section [5]) of the successive 
terms of the expansion with the expected expansion (see Section [6} near the critical point, 
one gets the periodic functions B\ and Bi- 

Starting from the map (jl.ip and defining the successive terms of the expansion of F(*)( x) 
in powers of e by 



F^\x) = x pn exp(ep n Y n (x) + e 2 p n Z n (x) + 0(e 3 )) . 
Then by replacing in the expression (jl.ip one gets 

Y n (x) = Y n _ l (x)+p- n G{x pn ' 1 ) , 

Z n (x) = Z n -x(x) - ^-G 2 (V"- 1 ) +-G' (x^ 1 ) x^ 1 Y n . x {x) . 



(4.1) 



(4.2) 
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One can show by recursion that 

Y n (x) = l -Y,^ G K) = m - ^J. , (4.3) 

fc=0 

where we recall U (x) from (|3,6p 

1 00 

U{x) = -Y,P~ k G(xV k ) , (4.4) 

and this leads to 



^ k=0 



i> n (x) = \ogx + eY n {x) + e 2 Z n (x) + 0(e 3 ^ 



n-l 



fc=o (4.5) 



_fc 



2 n ~ l / {7 X P 

For large n, assuming G(-) is such that the sums converge, one gets 

^oo (x) = logx + e C/(x) + e 2 (-^^ + V 2 (x) - Vg(a 

where 



1 oo 

y k=0 
-. oo 

y 2 (x) = -J2 xpk G ' { xpk ) > (4-7) 

^ fc=0 
1 oo 

y 3 (x) = ~Y,P~ k x pk G'(x pk ) u(x p 



Since G(-) is bounded, the first series converges. Moreover if the series for V^(-) converges, 
so does the series defining V^-). Therefore we require the converges of the series defining 
and this amounts to an additional assumption on the large x behavior of G"(x): for 
example, it suffices that |G"(x)| < l/x(logx) a ) for some a > 1 and x large. 

5. Singular behavior of the sums (|3.6p -(|4.7 p 

In this section we try to extract the singular behavior, as x \ 1, of the sums U (x), Vi(x), V^(x) 
and V 3 (x) defined in $3ft - $ZJ) . 

For x sufficiently close to 1 one can write 

log (logx) = -(n + y)logp, (5.1) 

np is an integer where < y < 1. Then as x \ 1, the integer hq — >■ 00. 

We will also assume that as in (ll.9p . one has G(l) = 0, and since G(-) is smooth 

G{x) = G'{l)(x - 1) + G"(l)^^ + ... (5.2) 
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5.1. Analysis of U(x) as x \ 1. One can easily see, using G(l) = and the definition 
(|5.ip of hq, that U{x) defined in (I3.6P can be written as 



U(x) = ^ ( J2 P y - k G (V^) + £ [ p y~ k G - G'(l 

^k=0 k=— no 

logx ( log ( log x) 



V 1o sp / 



p 

Then, by using (|5.2p . we see that tlq — >• oo as x \ 1 and 

C/(x) - 17* (x) = O (log 2 x) , (5.4) 

where U*(x) is defined as 

[/* (x) = - (log x) (log(log x)) + B 1 (y) log x , (5.5) 

where -Bi(-) is the periodic function {B\{y) = B\(y + 1)) defined in (|3.4I) . 

A few remarks: 

(1) One can show that 

J7(x) = -U(x p ) + -G(x) , (5.6) 
V P 

U*(x) = -U*(x p ) + -G'(l)\ogx. (5.7) 
V P 

(2) One can also evaluate the correction U(x) — U*(x) and get for x \ 1 

U(x) = U*(x) - G "M ± log 2 (x) + O (log 3 (x)) . (5.8) 

(3) One can check directly from the full expression (|3,4h of B\(y) that 

B^y + l) = B 1 (y). (5.9) 



5.2. Analysis of V±(x) as x \ 1. Using the definition f)5. 1 1) of tiq and the behavior (|5.2p 
of Cr(x) when x \ 1 one can rewrite the sum in the first line of (|4.7p as 

, oo 

V&Xz) = ^ £ p^ fc G 2 (e^*) , (5.10) 

^ k=~n 

and as n — )• oo in the limit x \ 1 one gets 

Vi(x) ~Ci(y)logx , (5.11) 
where Ci(-) is the periodic function (Ci(y + 1) = Ci(y)) defined by 

Ci(y) = \ £ ( eP " y ) • ( 5 - 12 ) 

fc=— oo 
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5.3. Analysis of V2(x) as x \, 1. Using again the relation (|5.1|) between x and no one 
can write the sum in the second expression in (14 .7h 



V 2 (x) 



+ ]T ( epk ~ V ) ~ C'(l) 



k=— no 



G'(l) 
V 



oo —i 



.k=0 



k=—oo 



log(logx) G"(l) 

-G (1) —y. 



plogp 



V 



One then finds 



V 2 (x) = -G'(l) l0g( ?° gx) + C 2 {y) + O(logx) , 



plogp 

with with the periodic function C 2 {-) given by 



C 2 (y) = - 
P 



oo —1 

. yG ' { l ) + Y j eP k - v G'(e^~ v )+ Y [e pk ~ V G'(f k - y )-G\l) 

k=0 k=—oo 



(5.13) 
(5.14) 

(5.15) 



Remark 5.1. One can easily see from the definitions (|3.6p and (|4.7p of U(x) and V 2 {x) 
that 

V 2 (x) = xU'{x). (5.16) 

This implies that 

B'M G'(l) 



C 2 {y) = B 1 {y) 



logp p logp 



(5.17) 



which can be also checked directly. 



5.4. Analysis of V 3 {x) as x \ 1. From the definitions in the third expression in (14. 7[) 
and (|5.1|) of the sum V3 and of no one can write 



V 3 (x) 



log X 

p 



(5.18) 



fe>0 



which can be rewritten using again (|5.1 
log x 



V 3 (x) 



P 



Y^ p ~ k +v e P V G' (e p v ) U (e p 



o k ~y 



k>0 



+ £ \ P~ k+y e pk ~ V G' (e? k - y ) U (e^) + ^^(A 



fc>— no 

it ft N \2 



P 



+ GW f -log(logs) \ / -logflogs) y | 1 



+ 



2p \ log p 

(yG'(\f 



logp 



V P 



(5.19) 
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Now from (|5.5|) one has for large negative k 

U (e^) = p k ~y (-^p-(k -y) + Bi(y)) , (5.20) 
so that one can take the no — > oo limit in the last expression of V%{x) to get 

V^) = - '° gXl f l0gX) B l{ y)Gi(l) 
piogp 

+ GW logx log(logx) fjogCggx)_A 
2p z logp \ logp / 

where the periodic function C^{y) is defined by 



C 3 (y) = - 



-^(^G'Cl) - „(„ + 1)^- + J2P~ k+y e pk - y G' (e^) U (e^ v 

P k>0 



+ £ L- fc +^ fc - B G / (e^)^(e^)+^^(A;-y)- J B 1 (y)G'(l)} 

fc=-oo ^ P J 

Here again one can check directly in (|5.22p that C^(y + 1) = C%(y). 



(5.22) 



6. Expected expression of Voo(-) 

Since here G(l) = and X* = 1, the expected critical behavior (|3.2p of ^oo(') should be 
of the form 

^o,(x) a (logxTB (- ^g^ ) +o((log*r +1 ) . (6.1) 
For e small one has 

F , (l)=p + e/ 1 , a = l + ea 1 + e 2 a 2 + ... and B(y) = 1 + eB^y) + e 2 B 2 {y) + . . . (6.2) 
with 

fi = G (1) , ai = : and a 2 = 7 r ? — , (6.3) 

p logp z(plogpj z 



and one gets from (16. ip 

^oo(x) = logx+ [Si (y) log x + ai log x log(logx))]e 

+ [B 2 (2/)logx + a 2 log(x) log(log(x)) + a 1 B 1 (y)logx log(log(x)) 



' e 2 + . . . (6.4) 



a , n x l2 , , iogx log(logx) ; 

+— logx log(logx) +/i B^y) 

2 p(logp) z 



where 

log (log x) 



\ogp 



(6.5) 



Comparing with (|4.6p and (16.4p with U(x), V\{x), V^x) and V^x) replaced by their 
estimates (foTol) . (jSTTTT) . (f5TT4|) and (f572T]) one then finds 

B 2 {y) = + Bi(y)C 2 (y) - C 3 (y) , (6.6) 
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which leads to (j3.5j) once the periodic functions C\,C2,C^ have been replaced by the 
expressions (f5TT2]l , (ETT5] ) and (1Q21) . 

7. Large p limit 

We consider now the the map (jl.lip : recall that q £ (0, 1) is fixed, as well as the 
polynomial P(-). When p is large x* = 1, F'(l) = pq + (1 — g)P'(l) ~ pq and one has 
(compare with (|1.7p ) 



a 



logp 



logg 



1 + 2 



logF'(l) logp 

For x close to 1, let us define z by 

i + (F'(i)) n{x) (x-r 

where n(x) is the integer such that 

p-* < z < F'{l)p-^ P ~°° qp 1 ^ 

and rj £ (0, 1) is some fixed number (for example rj = 1/3). 

The main idea is that for x close to 1, the map (jl.lip can be replaced by F(x) ~ 
1 + F'(l)(x — 1) for the first n(x) iterations, and by F(x) ~ qx p beyond the ra(x) + 2 nd 
iteration. Only the n(x) + 1 st iterate has to be calculated with the full expression (jl.lip . 

Therefore recalling that P(l) = 1 we have 



(7.1) 
(7.2) 

(7.3) 



pn{x)+m+l 



(X) ~ (< 



e 2 + l 



p"'-i 
i p-i 



and 



log F n ^ +m+1 (x) 

pn(x)+l+m 



-P 



-n(x) 



P 



log (q e z + 1 - q) + 



p — 1 p m (p — 1) 



log<7 



and by taking the m — > oo limit and using (F'(l)) a = p we obtain 



(x — l) a /p 
p V z 



log (q e z + l-q) + 



p-1 



Taking the large p limit we therefore obtain 



ipoo(x) ^ (x - iy 



log (g e z + 1 



One can then write 



that is 



log(x — 1) 



logF'fl) 



n(x) + y with rj < y < rj + 1 , 



n(x) 



log(x — 1) 
' loeF'(l) 



where [a] := max{n £ Z : n < o}, or equivalently (see (j7.2 

£ = (F'(l)P, 



(7.4) 
(7.5) 

(7.6) 

(7.7) 

(7.8) 
(7.9) 
(7.10) 



and note that in the large p limit z / oo (respectively z \, 0) for y < 1 (respectively 
y > 1). Then one finds, for fixed y / 1 



(x — l) a g y for ?7 < y < 1 , 
(x - l) ^ 1 "^ for l<y<l + r?, 



(7.11) 
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as claimed in (|3.14|) . 

At y = 1 the amplitude becomes discontinuous. One can resolve this discontinuity by 
analyzing, in the large p limit, the range y — 1 ~ (logp) -1 and one gets (|3.16|) from (|7.1[) 
and ([ZZD - 
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